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Abstract 

We study classical Hamiltonian systems in which the intrinsic proper time evolution pa- 
rameter is related through a probability distribution to the physical time, which is as- 
sumed to be discrete. In this way, a physical clock with discrete states is introduced, 
which presently is still treated as decoupled from the system. This is motivated by the re- 
cent discussion of "timeless" reparametrization invariant models, where discrete physical 
time has been constructed based on quasi-local observables. Employing the path-integral 
formulation of classical mechanics developed by Gozzi et al., we show that these determin- 
istic classical systems can be naturally described like unitary quantum mechanical models. 
We derive the emergent quantum Hamiltonian in terms of the underlying classical one. 
Such Hamiltonians typically need a regularization - here performed by discretization - 
in order to arrive at models with a stable groundstate in the continuum limit. This is 
demonstrated in two examples, recovering and generalizing a model advanced by 't Hooft. 

PACS: 03.65.Ta, 04.20.-q, 05.20.-y 



1 Introduction 

Since its very beginnings, there has been a long series of speculations on the possibility of deriv- 
ing quantum theory from more fundamental dynamical structures, possibly deterministic ones. 
Famous is the discussion by Einstein, Podolsky and Rosen. This lead to the EPR paradox, 
which in turn was interpreted by its authors as indicating the need for a more complete fun- 
damental theory pQ. However, just as numerous have been attempts to prove no-go theorems 
prohibiting exactly such "fundamentalism" . This culminated in the studies of Bell, leading to 
the Bell inequalities |2j. The paradox as well as the inequalities have come under experimen- 
tal scrutiny in recent years, confirming the predictions of quantum mechanics in laboratory 
experiments on scales very large compared to the Planck scale. 



* E-mail address: thomas@if.ufrj.br 



2 



However, to this day, the feasible experiments cannot rule out the possibility that quantum 
mechanics emerges as an effective theory only on sufficiently large scales and can indeed be 
based on more fundamental models. 

Motivated by the unreconciled clash between general relativity and quantum theory, 't Hooft 
has strongly argued in favour of model building in this context j3] (see also further references 
therein). He has shown in individual examples the emergence of the usual Hilbert space struc- 
ture and unitary evolution in deterministic classical models in an appropriate large-scale limit. 
Particular emphasis has been placed on the observation that while it is relatively easy to ar- 
rive at a Hilbert space formulation for classical dynamics, it is difficult to obtain emergent 
Hamiltonians with a spectrum bounded from below, i.e., having a well-defined groundstate. 

Various further arguments for deterministically induced quantum features have recently been 
proposed, for example, in Refs. jHEUHllZj, in the context of statistical systems, of considerations 
related to quantum gravity, and of matrix models, respectively. 

Our aim here is to contribute to this line of research by reporting a rather large class 
of classical deterministic systems which yield to a quantum mechanical description. This is 
based on our recent work on time-reparametrization invariant models [HUH!- In particular we 
have introduced the construction of a discrete physical time for such "timeless" systems, which 
will be the starting point of the present developments. In the following we briefly summarize 
essential aspects. 

This paper does not claim to have derived quantum mechanics from an underlying classical 
dynamics yet. However, we develop a formalism which may prove useful for this purpose, besides 
having potentially other interesting applications in statistical mechanics. We do illustrate our 
approach with first examples indicating the possible emergence of quantum mechanics from 
deterministic models. 1 

1.1 Discrete physical time in "timeless" classical models 

Similarly to common gauge theories, such as those of the standard model of particle physics, 
reparametrization invariant systems show invariance under a kind of gauge transformations. In 
the most general case of diffeomorphism invariant theories, such as general relativity or string 
theory, this amounts to invariance under general coordinate transformations. 

We limit ourselves to time-reparametrization invariance here, which for our purposes can 
be expressed as invariance of the dynamics under arbitrary transformations: 2 



Details of the corresponding constrained Lagrangian dynamics can be found in Refs. jSl E| for 
the respective models. 

Similarly as Gauss' law in electrodynamics, for example, a most important consequence of 
time-reparametrization invariance is a (weak) constraint, which states that the Hamiltonian has 
to vanish (on the solutions of the equations of motion). 3 Since the Hamiltonian commonly is 

1 A summary of this work was presented as invited talk at the International Workshop Trends and Perspectives 
on Extensive and N on- Extensive Statistical Mechanics in honor of Constantino Tsallis' 60th birthday, Angra 
dos Reis (Brazil), Nov.19-21, 2003 QUI. 

2 Some restrictions are imposed in Refs. jHUHj on physical grounds, such as differentiability and monotonicity 
of the function /. 

3 In the case of the free relativistic particle the invariance amounts to invariance under reparametrization of 
the proper time parameter and the ensuing constraint is the familiar mass-shell constraint on its momentum. 
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the generator of time evolution of the system, this is what has led to name this type of systems 
"timeless". A Newtonian external time parameter does not exist. This becomes problematic 
when trying to quantize such systems, since a standard Schrodinger equation does not exist ei- 
ther. In particular, the Wheeler-DeWitt equation, H\ip) = 0, epitomizes the intrinsic problems 
of quantum gravity, seen from this perspective. 

Numerous approaches have been tried to resolve this (in)famous "problem of time". For 
our purposes, it may suffice to mention one, in order to contrast it with our proposal. 

In Refs. [TTJ H2 EH] it has been assumed that global features of a suitably parametrized 
trajectory of the system, are accessible to the observer. This makes it possible, in principle, to 
express the evolution of an arbitrarily selected degree of freedom relationally in terms of others. 
Naively speaking, the question "What time is it?" is replaced by "What is the value x of ob- 
servable X, when observable Y has value y?" . Thereby the Hamiltonian and possibly additional 
constraints have been eliminated in favour of Rovelli's "evolving constants of motion". 

In distinction, we insist on a local description. We have shown that for a particle with 
time-reparametrization invariant dynamics, be it relativistic or nonrelativistic, one can define 
quasi-local observables which characterize the evolution in a gauge invariant way (HI E] • 

Essentially, we employ some of the degrees of freedom of the system to trigger a localized 
"detector" . This can be defined in an invariant way. It amounts to attributing to an observer 
the capability to count discrete events. - In passing we remark that a similar approach, avoiding 
the notion of time altogether and replacing it by counts of idealized coincidence detectors in 
phase space has recently also been put forth by Halliwell and Thorwart [H] . - In any case, the 
detector counts present an observable measure of time. This physical time is discrete. 

This result can be understood in a different way, by noting that our construction is practi- 
cally based on a Poincare section, or subsection thereof, which reflects an ergodic if not periodic 
aspect of the dynamics - quite analogous to a pendulum which triggers a coincidence counter 
each "time" it passes through its equilibrium position. Reparametrization invariance strongly 
limits the information which can be extracted from it with respect to a complete trajectory. 
This is the underlying reason that a physical time based on local observations (clock readings) 
necessarily is discrete. 4 

Independently of our physical motivation of discrete time, we remark that the possibility 
of a fundamentally discrete time (and possibly other discrete coordinates) has been explored 
before, ranging from an early realization of Lorentz symmetry in such a case to detailed 
explorations of its consequences and consistency in classical mechanics, quantum field theory, 
and general relativity ^21 El EH] ■ Further recent developments are discussed in the review of 
Ref. ^H]; particularly in relation to fundamental or induced violations of Lorentz symmetry, 
which are believed to come within experimental reach in the near future. 

So far, however, no classical physical models giving rise to such discreteness were proposed. 
Quantization as an additional step - which results in discreteness of coordinates in various cases 
- has always been performed as usual. 



4 This does not conflict with certain cosmological models where a scalar matter field coupled to gravity is 
invoked as a continuous time variable. Here, the need to know the scalar field globally can be circumvented by 
assuming its homogeneity. 
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1.2 Does discrete time induce quantum mechanical features? 

It is the purpose of the present study to reach a qualified 'Yes', answering this question. We 
hope this will contribute to the study of potentially deterministic substructures leading to 
quantum mechanics as an emergent theory. 

Based on the findings of Refs. [HUH], respectively, we have argued there that those discrete- 
time models can be mapped on a cellular automaton studied by 't Hooft before j3] . With the 
help of the algebra of SU(2) generators, it has been shown that these models actually reproduce 
the quantum mechanical harmonic oscillator in a suitably defined large-scale limit. 

We will come back to variants of the earlier results in Section 6. However, in the following 
main parts of this paper, we attempt to show more generally that due to inaccessability of 
globally complete information on trajectories of the system, the evolution of remaining degrees 
of freedom appears as in a quantum mechanical model when described in relation to the discrete 
physical time. 

We may call this "stroboscopic" quantization: when a continous physical time is not avail- 
able but a discrete one is - like reading an analog clock under a stroboscopic light - then 
states of the system which fall in between subsequent clock "ticks" cannot be resolved. (Of 
course, evolution in the unphysical parameter time is continous in the constrained Lagrangian 
models we refer to.) Such unresolved states form equivalence classes which can be identified 
with primordial Hilbert space states [3 HJ |H1 El • The residual dynamics then has to describe the 
evolution of these states through discrete steps. Under favourable circumstances, this results 
in unitary quantum mechanical evolution, as we shall see. 

Our present aim is to show that this occurs quite generally in classical Hamiltonian systems, 
if time is discrete. We will presently simplify the situation by assuming that the physical time 
can be related by a probability distribution to the proper time of the equations of motion. 
Explicit examples of such behaviour can be found in Refs. [HI E|, when the clock degrees of 
freedom evolve independently of the rest of the system, apart from the Hamiltonian constraint. 

Thus, while the investigation of the coupled system-clock dynamics is presently under study 
[20J, here we make the approximation that corresponding "backreaction" effects are small. This 
feature, besides characterizing a good clock, may serve as an appropriate simplification in our 
exploratory steps. In the concluding section, we will briefly comment about extensions, where 
the prescribed probabilistic mapping of physical onto proper time shall be abandoned in favour 
of a selfconsistent treatment. A closed system, of course, has to include its own clock, if it is 
not entirely static, reflecting the experience of an observer in the universe. 

Finally, in order to put our approach into perspective, we remark that there is clearly no need 
to follow such construction leading to a discrete physical time in ordinary mechanical systems 
or field theories, where time is an external classical parameter. However, assuming for the time 
being that truly fundamental theories will turn out to be diffeomorphism invariant, adding 
further the requirement of the observables to be quasi-local, 5 when describing the evolution, 
then such an approach seems natural, which promises to lead to quantum mechanics as an 
emergent description or effective theory on the way. 

The paper is organized as follows. Section 2 provides a selfcontained brief summary of the 
path-integral formulation of classical mechanics. We will employ this as a convenient tool 
to formulate our approach. In later developments one might also introduce a cooresponding 



5 For example, a closed string loop representing a fundamental length might define the ultimate resolution of 
distance measurements and, thus, condition the notion of locality. 
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operator formalism, for example, considering as a starting point the early work of Koopman 
and von Neumann [21] . 

In Section 3 the relevant to-be-quantum states are introduced, i.e. equivalence classes of 
states of the underlying classical system, as we mentioned. Their evolution is studied in Sec- 
tion 4, leading to an emergent Hamilton operator and associated Schrodinger equation under 
circumstances to be discussed there. In Section 5 the calculation of physical time dependent 
observables, or expectation values of corresponding operators, is related to the states. 

Section 6 presents some simple examples of the emergent Hamiltonians and the calculation of 
their spectra. We show that - under the present simplifying assumptions - these operators need 
to be regularized, in order to represent acceptable quantum models, with a stable groundstate in 
particular. Here we achieve this by discretization. We speculate that its apparent arbitrariness 
might be removed by a future selfconsistent treatment of clock-system interactions, which will 
lead to dissipative effects on small scales and could define a unique quantum system for each 
classical model with discrete time. The concluding Section 7 presents a brief summary of the 
presented work and points out some open problems. 



2 Classical mechanics via path-integrals 

Classical mechanics can be cast into path-integral form, as originally developed by Gozzi, Reuter 
and Thacker [22], arid with recent addenda reported in Ref. [22] • While the original motivation 
has been to provide a better understanding of geometrical aspects of quantization, we presently 
use it as a convenient tool. We refer the interested reader to the cited references for details, 
on the originally resulting extended (BRST type) symmetry in particular. Here we suitably 
incorporate time-reparametrization invariance, assuming equations of motion written in terms 
of proper time (as in our earlier examples [HIE])- 

Let us begin with a (2n)-dimensional classical phase space Ai with coordinates denoted 
collectively by ip a = (g 1 , . . . , g n ; p 1 , . . . ,p n ), a = 1, . . . , 2n, where q,p stand for the usual 
coordinates and conjugate momenta. Given the proper-time independent Hamiltonian H(ip), 
the equations of motion are: 

|y , (2) 

where uj ab is the standard symplectic matrix and r denotes the proper time; summation over 
indices appearing twice is understood. 

To the equation of motion we add the (weak) Hamiltonian constraint, Ch = H(tp) — e ~ 0, 
with e a suitably chosen parameter. This constraint has to be satisfied by the solutions of the 
equations of motion. Generally, it arises in reparametrization invariant models, similarly as 
the mass-shell constraint in the case of the relativistic particle 0|. It is necessary when the 
Lagrangian time parameter is replaced by the proper time in the equations of motion. In this 
way, an arbitrary "lapse function" is eliminated, which otherwise acts as a Lagrange multiplier 
for this constraint. 

We remark that field theories can be treated analogously, considering indices a, b, etc. as 
continuous variables. 

Starting point for our following considerations is the classical generating functional, 

Z[J] = j^Dtp 5[<p a (r) - exp(* J dr J a <p a ) , (3) 
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where J = { J a =i,...,2n} is an arbitrary external source, S[-] denotes a Dirac ^-functional, and ipd 
stands for a solution of the classical equations of motion satisfying the Hamiltonian constraint; 
its presence is indicated by the subscript "if" on the functional integral. It is important to 
realize that Z[0], as it stands, gives weight 1 to any classical path satisfying the constraint and 
zero otherwise, integrating over all initial conditions. 

Using the functional equivalent of 5(f(x)) = \df/dx\~* ■ 8(x — xq) , the ^-functional under 
the integral for Z can be replaced according to: 

5[ V %t) - ^(r)] => 5[d T tp a - uj ab d b H] det[6 a b d T - u ac d c d b H] , (4) 

slightly simplifying the notation, e.g. d b = d/d(p b . Here the modulus of the functional deter- 
minant has been dropped (221 122] ■ 

Finally, the 5-functionals and determinant are exponentiated, using the functional Fourier 
representation and ghost variables, respectively. Thus, we obtain the generating functional in 
the convenient form: 

Z[J] = J^VtpVWcVc exp (i J dr(L + J a (p a )) , (5) 

which we abbreviate as Z[J] = f H T>§ exp(i J drLj). The enlarged phase space is (8n)- 
dimensional, consisting of points described by the coordinates ((p a , X a , c a , c a ). The effective 
Lagrangian is now given by j221 |2*H] : 

L = X a (d rV a - u ab d b H) + ic a (5 a h d T - u ac d c O b H) c b , (6) 

where c a , c a are anticommuting Grassmann variables. We remark that an entirely bosonic 
version of the path-integral exists [23J . 

This completes our brief review of how to put (reparametrization invariant) classical me- 
chanics into path-integral form. 

3 From discrete time to "states" 

We recall from the motivation provided in Section 1.1 that discrete physical time t is constructed 
and based on the counting of suitably defined incidents. In particular, we have coincidences in 
mind, "when" points belonging to the trajectory of the system coincide with the position of 
an idealized detector. For concrete realizations of this procedure and study of its invariance as 
well as further intrinsic properties we refer the interested reader to Refs. [HI EI , and to Ref. (Ej 
for a similar construction. 

Thus, physical time is measured by a nonnegative integer multiple of some unit time, t = nT, 
and is read off from a sort of localized track counting device. 6 For our present study details 
of the clock construction are of secondary importance. Rather we investigate the consequences 
of the discreteness of physical time. In particular, we want to demonstrate the existence of 
quantum mechanical features which may derive from it. 

In any case, then, the proper time r parametrizing the evolution should be calculable in 
terms of t. As can be observed, for example, in our earlier numerical simulations, this can be 

6 In our earlier detailed examples we were motivated by the at present actively researched separation in higher 
dimensional (cosmological) models of bulk and brane degrees of freedom: interactions between both kinds of 
matter could lead to a more detailed picture of such countable incidents providing the basis of physical time. 
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a formidable task, depending on the dynamics governing the underlying classical system [SJ E] • 
Therefore, an analytic approach for interacting clock-coupled-to-system models is beyond the 
scope of this work and considered elsewhere [2Uj . 

We will assume that the backreaction and memory effects which generally result from the 
coupling between clock and system are small and negligible for our present purposes. Then the 
relation between the discrete physical time t and the proper time r of the equations of motion, 
Eq. (0), can be represented by a time independent normalized probability distribution P: 

P(r;t) = P(r-t) = exp(-S(r-t)) , J dr P(r; t) = 1 . (7) 

Note, in particular, that the perfect clock described in this way does not age with physical 
time. Eventually, we will also invoke the limiting case, P(r;t) = 6(r — t), i.e. a deterministic 
mapping between t and r. 

We remark that in the present situation, the Hamiltonian constraint needs to be applied 
to the system degrees of freedom only, while generally clock plus system are constrained as a 
whole. This is exemplified in all detail in our toy models [HE]- 

Correspondingly, we introduce the modified generating functional: 

Z[J] = f dndy [ V<$> exp (i f Tf dr Lj - S(n - U) - S(r f - t f j) , (8) 

instead of Eq. (j5J), and using the condensed notation introduced there. In the present case, Z[0] 
sums over all classical paths satisfying the constraint with weight Pfa] ti) -P{Tf, tf), depending 
on their initial and final proper times, while all other paths get weight zero. In this way, the 
distributions of proper time values Tij associated with the initial and final physical times, U 
and tf, respectively, are incorporated. 

Next, we insert 1 = / drP(r; t) into the expression for Z, with an arbitrarily chosen physical 
time ti < t < tf. This leads us to factorize the path- integral into two connected ones: 

S(Tf - t f )) 

S(r, - U)) 

■ n -*><(*-)) > ( 9 ) 

a 

where "<" and ">" refer to earlier and later than r, respectively. 

The ordinary ^-functions assure continuity of the classical paths in terms of the coordinates 
q a , a = 1, ... ,n and momenta p a , a = 1, . . . ,n, at the (distributed) proper time r. This is 
necessary, in order to avoid a double-counting as compared to the original expression for Z, 
Eq. (jSJ). Since any given classical path beginning at ti and ending at tf (with associated values 
Tij) is cut in two parts, without the continuity condition, the continuing part (">") would again 
have arbitrary initial conditions at the instant r, unlike the corresponding original path where 
this is excluded by the Hamiltonian flow. To put it differently, without continuity, originally 
N independent paths would be factorized erroneously into N 2 independ ones. 

Another remark is in order here. Due to the presence of the Hamiltonian constraints on 
both parts of the cut trajectory, one of the (^-functions is redundant. The resulting 6(0) may 
eventually be absorbed into the normalization of the states, which will be introduced now. 



Z[J] = j dr P(r;t) ■ J dr f J exp (i 1 dr' Lj — 

■ JdT t J^V^ < exp (i £ dr" Lj 
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Exponentiating the (^-functions via Fourier transformation, the generating functional can 
indeed be interpreted as a scalar product of a state and its adjoint: 

Z[J] = J drdic P(T)(t;t f \T,Tr)(T,Tr\t\U) , (10) 
in "t, 7r-representation" ; here d7r = Yl^diia /2ti). The state is defined by the path-integral: 

(r,n a \t;t Q ) = Jdr' J h V1> exp(i£ +t dr"Lj-S{r' -t )+i7i aV a {T + t)) , (11) 
and, similarly, the adjoint state: 

(t;t \r,7c a ) = J dr' J h V$ exp (i £ dr"Lj - S(r' - t ) - m a (p a {r + tj) , (12) 

where to represents an arbitrary physical time with which the respective paths begin, in the 
case of states, and end, in the case of adjoint states. 

We observe a redundancy in designating the states, which depend on the sum of proper and 
physical time only. This arises by the shift of r leading to Eq. (jlUJl . noting that the probability 
distribution P is not explicitly depending on the physical time, as we discussed. 

Furthermore, there is a symmetry between states and adjoint states: 

(t;to|r,7r) = (r,7r|t;to>* , (13) 

if referred to the same reference time to, which is familiar in Hilbert space. 
Finally, the scalar product of any two states can be defined by: 

(t 2] t f \tuU) =Z[J]- 1 |drd7rP(r)(t 2 ;t / |r,7r)(r,7r|t 1 ;t i ) . (14) 

Of course, we have {t;tf\t;tj) = 1, corresponding to Eq. (J0J), while for t 2 ^ t\, the states 
generally are not orthogonal. 

Closing this section, we point out that here the path-integrals have always been considered to 
include all classical paths. However, they can be restricted by imposing (final) initial conditions 
in the case of (adjoint) states. Since this has no effect on the following study of the time 
evolution of a generic state, there is presently no need to make explicitly use of this. 

4 Unitary evolution 

In order to learn about the time evolution of generic states (in the absence of a source, J = 0), 
we proceed similarly as in Eqs. (J5j)- (fTTj) . 

Suitably inserting "1" , as before, and splitting the path-integral (at time t) which defines a 
state (at time t'), we obtain: 

( T ',n'\t') = J drdir P{T)(T',ir'\U{t',t)\T,7c)(T,7c\t) , (15) 

where henceforth we suppress to indicate the reference time to, see Eq. (jllj) . in order to simplify 
the notation. Written as a matrix element of an evolution operator U(t',t), the kernel which 
appears is: 

r rr'+t' 

(T',7['\U(t',t)\T,7i) = J h 'D§ exp (i J dr"L + m ■ ip{r' + t') - m ■ ip{r + t)) , (16) 
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with the functional integral over all paths running between r + t and r' + t', subject to the 
Hamiltonian constraint; we abbreviate: 7r • (p = n a (p a . 

Then, first of all, it is straightforward to establish the following composition rule: 

U(t",t').U(t',t) = U(t",t) , (17) 

written here in the way of a matrix product which is to be interpreted as integration over 
intermediate variables, say dr', dn' , with appropriate weight factor P(t'). Then, the first 
integration, say over dr', effectively removes a "1" - such as we inserted before, in order to 
factorize path-integrals. The second, say over dn' reinstitutes (^-functions - such as those which 
first made their appearance in Eq. Q - which serve to link the endpoint coordinates of one 
classical path to the initial of another. This leads to the result of the right-hand side. 

Since the Hamiltonian constraint is a constant of motion, there is no need to constrain 
the path-integral representing the evolution operator. Integrating over intermediate variables 
removes all contributions violating the Hamiltonian constraint, provided we work with properly 
constrained states. This will be further discussed in the following section. 

We observe in Eq. (JT5|) that the physical-time dependence of the evolution operator (on t' 
and t) amounts to translations of proper time variables (r' + t' and T+t) in the matrix elements. 
This simplicity, of course, is related to the analogous property of the states, see Eq. (jllj) . which 
we discussed in the previous Section 3. 

Tracing the steps backwards which led from Eq. (jSJ) to Eq. (JHJ in Section 2, we can similarly 
rewrite here the functional integral of Eq. (|16|) : 

(r',n'\U(t',t)\r,n) = J V<p % a (f ) - (p a d (f )] exp (in' ■ <p(r' + t') - m ■ <p(r + t)) , (18) 

where the paths parametrized by f run between r + 1 and r' + 1' and all initial conditions are 
integrated over. 

Fixing momentarily the initial condition of the paths, <p(r + 1) — (p^ we select a particular 
path contributing in Eq. ()18j) . By integrating over all initial conditions in the end, / dipi, we 
recover the full expression. Again, this amounts to a suitable insertion of "1", a familiar 
procedure by now, which allows us to manipulate the above equation as follows: 

(r',7r'\ U(t', t)\r,7r) (19) 

= jv<pj difi 5(<p(r + t)- (pi) % a (f ) - ^(f)] exp (V ■ cp{r' + t') - iir ■ cp{r + £)) 

= J dtp,, exp (in' ■ <p f - in ■ ^) J Vtp 5[tp a (f) - (p a d (f)) , 

where cpf denotes the value of (p(f) at the endpoint r' + t' of the path singled out by the initial 
condition (p(r + t) = (pi. This value is determined by the classical equations of motion and 
will be calculated shortly. Since the remaining functional integral equals one - only the one 
path with fixed initial conditon and weight 1 is contributing - we first obtain the intermediate 
result: 

(r',iT'\U(t',t)\T, 7r) = J difi exp (ire' ■ <p f - in ■ ip^ , (20) 

where d<^ = Ulli d^. Of course, the time dependent relation between (p$ and (pi has to be 
enforced, in order to make this condensed expression explicit. This we do next. 
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We employ the Liouville operator, 

where is the symplectic matrix, stemming from the Hamiltonian equations of motion As 
is well known, the Liouville operator allows to generate a classical solution of these equations 
at any finite time, starting with a given initial condition. Applying this for our purposes yields: 

(p f = (p(r' + t') = exp[£(r' + t' - r - t)]cp(r +t)= exp[£(r' + t'-r- t)](p { . (22) 

Then, inserting into Eq. ()20|) . we obtain the simple but central result: 

(T',n'\U{t',t)\r,n) = J d<p exp (in' ■ (exp[£(r' + t' - r - t)](p) - in ■ cpj (23) 

= £(tt', tt; t' + t' - r - t) , (24) 

where the by now superfluous subscript "z" has been omitted. 

Using Eq. ()23|) . one readily confirms Eq. p7|) once again. In particular, then U(t\t') ■ 
U(t'\t) = U(t\t), which is not diagonal, in general, in this r, 7r-representation. We have: 
U(t', tt'; t\r, ?r; t) = £(n', n; t'-t). 

In order to proceed, we consider the time dependence of the evolution kernel £. Beginning 
with Eq. (j2U), one derives the equation: 

id T £(n', tt; t) = — J dcp exp (in' ■ <p(r) — n ■ y?) tt' • (d T (p(r)) 

= — J d(p exp (in' ■ (p(r) — n ■ (p^j n' ■ uo ■ —H((p(r)) 

= H{n\-id nl )£{n\n]T) , (25) 

with the effective Hamilton operator. 

— d 

H(tt, -idn) = -tt ■ uj ■ — Hicpjl^-i^ . (26) 

Here we also used Eq. (J2~2"j) in the first step, the equations of motion (J2J) in the second, and 
suitably pulled the factor following the exponential out of the integral at last. The initial 
condition, 

£( 7r ', n- 0) = {2n) 2n 5 2n (n' - n) , (27) 
can be read off from Eq. ()23j) . Integrating Eq. (|25jl. we immediately obtain: 

£(Tr', tt; t) = (2n) 2n expf-ir^Tr', -id^)]5 2n {n' - n) , (28) 

taking the initial condition into account. 

This result, which yields the evolution operator U by Eq. (|24j) . finally allows to describe the 
evolution of a generic time dependent state, \^/(t)), which takes place in one discrete physical 
time step (unit T). Using Eq. (fTK|) . we calculate: 

(T',Tr'\ty(t + T)} = J drdn P(t)£(tt',tt;t' + T - r)(r,7r|*(t)) 

= Jdr P(T)exp[-i(T' + T -T)H(n',-id n ,)}(T,n'\ty(t)} , (29) 
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where the integration over d7r = Y\ a {&K a /2n) has been carried out with the help of the 8- 
functions of Eq. ()28|) . This equation appears like a discrete time Schrddinger equation and 
presents our main result of this section. However, some qualifying remarks are due here. 

4.1 Discussion 

In order to facilitate the investigation of the properties of Eq. (j2T?j) . let us assume that a generic 
time dependent state shares the following (proper) time translation property with the state 
defined in Eq. (fTTj) : 

(t, tt|# (*)) = (r + t, 7r|*(0)) = (0, ir\^{t + r)) . (30) 

This is the case, for example, if \^(t)) is defined like \t) there, however, with some particular 
initial conditions for the paths contributing to the functional integral, or if it is a superposition 
of such states. 

Using this property and shifting and renaming variables, the Eq. ()29j) can be rewritten: 

(0,7r'|tf(r' + T)) = J dr P(t -t)exp[-i(r' + T - r)7?(7r', -^)](0,7r'|^(r)) . (31) 

Despite that a state needs to exist only on discrete values t of the physical time, the Eqs. (JBUj) 
or ()31|) require the corresponding "wave function" in r, ^-representation to be analytically 
continued to arbitrary real values of the time argument. 

Considering the deterministic limiting case, P(r — t) = 8(r — t) we obtain directly from 
Eq.flHH): 

(0, Tr'ltf (t' + T)) =exp[-z(r , + T-t)7?(7r',-z^)](0,7r'|^(t)) . (32) 

Thus, the usual formal solution in terms of an exponentiated Hamilton operator of a standard 
Schrddinger equation is recovered, i.e., the Eq. (J2~9"j) here is equivalent to such an equation. 

In more general cases, for example, with a Gaussian probability distribution, P(t — t) oc 
exp[— 7(t — t) 2 }, one might suspect that the evolution operator on the right-hand side of Eq. (J29|) 
or (|31|) is not unitary, even if 7i is hermitean. However, it is easily seen that the Ansatz, 

(r,7r|*(t)) =expH(r + 0^,-^)1(0, tt|*(0)) , (33) 

solves Eq. f)29j) also for any normalized distribution with P(t; t) = P(r — t). In this sense, the 
Eq. f)29j) is indeed formally equivalent to a Schrddinger equation. We will study this in more 
detail in Section 6, paying attention to the nonstandard form of the effective Hamiltonian (|26|) 
in several examples. 

The above Ansatz suggests to introduce stationary states defined by the following relations: 



(r,7r|* n (t)) = exp(-z^)(r,7r|* n (0)) = exp(-iS n (t + r))(0,7r|* n (0)) (34) 

= expH£„(t + T))(7r|^„> , (35) 

in accordance with Eqs. (JBPj) . and assuming: 

7?(7T, -1^)1(0, 7T|* n ) =E n (0,7l\^ n ) , (36) 



with a discrete spectrum {E n }, in order to be definite. 



12 



Finally, we remark that there is no % in our equations. If introduced, it would merely act 
as a conversion factor of units. On the other hand, there is the intrinsic scale of the clock's 
unit time interval T. The significance of this can be analyzed in a treatment where clock and 
system are part of one universe and interact [20 . 

Before we will illustrate in some examples the type of Hamiltonians and stationary "wave 
equations" that one obtains, we have to first address the classical observables and their place 
in the emergent quantum theory. In particular, we need to implement the classical Hamil- 
tonian constraint, which is an essential ingredient related to the gauge symmetry in a time- 
reparametrization invariant classical theory. 

5 Observables 

It follows from our introduction of states in Section 3, see particularly Eqs. (|8j)- (fT2j) . how the 
classical observables of the underlying mechanical system can be calculated. Considering ob- 
servables which are function(al)s of the phase space variables if, the definition of their expec- 
tation value at physical time t is obvious: 

= J dTP(T;*)0[-<^] log Z[ J] | J=0 (37) 

= Z^JdrdTT P(r-t)(0|r,7r)O[v9(r)](r,7r|0) (38) 

= Z' 1 J drdTT P(r)(t\T,7r)0[<f(T + t)](r,7r\t) (39) 

= Z^Jdrdn P(T)(t\T,7r)0[-idv](T,7r\t) (40) 

= (¥(f)|d[d|¥(f)> , (41) 

where Z = Z[0], all states refer to J = as well, and where: 

6[<p] = O[0\ , , (42) 

in r, 7r-representation. In Eqs. (|38|) - (|39j) the notation is symbolical, since the observable should 
be properly included in the functional integral defining the ket state, for example. Furthermore, 
in the last Eq. (|41|) . the preceding expression is rewritten for the case of a generic state 
(J = 0), with the scalar product to be evaluated as in (j4~Uj) . or as defined in Eq. flT3} before. 

Thus, a classical observable is represented by the corresponding function(al) of a suitably 
defined momentum operator. Furthermore, its expectation value at physical time t is repre- 
sented by the effective quantum mechanical expectation value of the corresponding operator 
with respect to the physical-time dependent state under consideration, which incorporates the 
weighted average over the proper times r, according to the distribution P. Not quite surpris- 
ingly, the evaluation of expectation values involves an integration over the whole r-parametrized 
"history" of the states. 

Furthermore, making use of the evolution operator U of Section 4, in order to refer observ- 
ables at different proper times to a common reference point r, one can construct correlation 
functions of observables as well, similarly as in Ref. [3], for example. 
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The most important observable for our present purposes is the classical Hamiltonian, H(tp), 
which enters the Hamiltonian constraint of a classical reparametrization invariant system. It 
is, by assumption, a constant of the classical motion. However, it is easy to see that also 
its quantum descendant, H((p) = H((p), is conserved, since it commutes with the effective 
Hamiltonian of Eq. (|26j) : 

— — d d 

[H,H\ = H(-id v ) tt-uj- — i/XvOI^-ifi* -TT-w- q- #(v0l<p=-iflw Hi-id*) 

d d 
= —H{ip)\ v= ^-u-—H{<p% = ^ = , (43) 

due to the antisymmetric character of the symplectic matrix. Therefore, it suffices to implement 
the Hamiltonian constraint at an arbitrary time. 

Then, the constraint of the form Ch = H(ip) — e ~ could be incorporated into the defi- 
nition of the states in Eq. (jll|) by including an extra factor 8{Cjj) into the functional integral, 
and analogously for the adjoint states. Exponentiating the (^-function, we can pull the expo- 
nential out of the functional integral, similarly as before. Thus, we find the following operator 
representing the constraint: 

C = J d\ exp (i\{H{<p) - e)) = S(C H ) , (44) 

which acts on states as a projection operator. It admits in the functional integral representing 
a state only those paths that obey the constraint; in particular, see Eq. ((lip, it enforces the 
constraint at the time r + 1. Similarly, a projector can be included into the definition of the 
generating functional, Eq. (|SJ), in order to represent the Hamiltonian constraint. 

Supplementing Eqs. (|37|) -()41 |) by the insertion of the Hamiltonian constraint, the properly 
constrained expection values of observables should be calculated according to: 

(0[<p]it) B ={^\d[<p]d\9{t)) , (45) 

which, in general, will deviate from the results of the previous definition. 

Finally, also the eigenvalue problem of stationary states, see Eqs. (|34| ) -([36 p . should be stud- 
ied in the projected subspace: 

nd\v) = e6\v) , (46) 

to which we shall return in the following section. 



6 Examples of quantum systems with underlying 
deterministic dynamics 

The purpose of the following examples is to illustrate how the quantum mechanics works in 
practice which emerges from various deterministic classical systems along the lines presented 
in Sections 2-5 before. In particular, we solve the stationary Eq. (|36jl with the effective Hamil- 
tonian H of Eq. pBjl . 

We shall see, however, that acceptable quantum models with a stable groundstate can only 
be arrived at in this way, if a regularization of the respective Hamiltonian and subsequently 
a suitable continuum limit are performed. In a sense, the spectrum of the Hamiltonian (|26|) 
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is too rich, it admits additional unphysical states, which have to be eliminated. While this 
approach covers a large number of models, the meaning of and possible restrictions on the ad 
hoc adopted regularization certainly deserve further study 

Our first example, starting with a "timeless" classical harmonic oscillator, mainly serves to 
demonstrate that the present general formalism allows to recover results of 'tHooft's cellular 
automaton model [Sj. The second model, employing the parameterized classical relativistic 
particle, is considered by many as a testing ground for any techniques developed to deal with 
reparametrization invariant theories in general, like general relativity or string theory. It leads to 
an interacting quantum model, provided we judiciously choose the arbitrary phases introduced 
by the regularization. 

6.1 Quantum system with classical harmonic oscillator beneath 

In principle, all integrable models can be presented as collections of harmonic oscillators. There- 
fore, we begin with the harmonic oscillator of unit mass and of frequency Q. The action is: 

S = /dt(^-(^) 2 -^(fiV-2e)) , (47) 

where A denotes the arbitrary lapse function, i.e. Lagrange multiplier for the Hamiltonian 
constraint, and e > is the parameter fixing the energy presented by this constraint. 

Introducing the proper time, r = / dt A, the Hamiltonian equations of motion and Hamil- 
tonian constraint for the oscillator are: 

8 T q = p , d T p = -Q 2 q , (48) 
^ + n 2 q 2 )-e = , (49) 

respectively. 

Comparing the general structure of the equations of motion, Eq. (J2J), with the ones obtained 
here, we identify the effective Hamilton operator, Eq. (|26jl. while the constraint operator follows 
from Eq. (jjjj) : 

H = -{n q {p p - fi 2 vr p ^) = -ir q (-id^) + Q 2 Tr p (-id^ q ) , (50) 
C = 6(0* + tf(p 2 - 2c) = 5{dl + tfd 2 q + 2c) , (51) 

respectively. Here we employ the convenient notation ip a = (ip q ;(p p ), and correspondingly 
7r a = (n q ;7i p ), d% = (d nq ;d Wp ). Further simplifying this with the help of polar coordinates, 
iT q = —Qp cos and tt p = p sin 0, we obtain: 

n = nl z = -ittd^ , (52) 

C = 5(A 2 + 2e) = 5(d p 2 + p~ x d p + p- 2 d} + 2e) , (53) 

where L z denotes the z-component of the usual angular momentum operator and A 2 the Lapla- 
cian in two dimensions. 

We observe that the eigenfunctions of the eigenvalue problem posed here factorize into a 
radial and an angular part. The radial eigenfunction, a cylinder function, is important for 
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the calculation of expectation values of certain operators and the overall normalization of the 
resulting wave functions. However, it does not influence the most interesting spectrum of the 
Hamiltonian. 

In order to proceed, we discretize the angular derivative, Eq. (J52J) . In the absence of the full 
angular momentum algebra, we would otherwise encounter a discrete yet unbound spectrum, 
lacking a groundstate. 

While we will mostly choose to work with an asymmetric discretization ( Case A ), we will 
here also show the symmetric discretization ( Case B ), in order to appreciate the differences, if 
any. In any case, the spectrum should and will turn out to be independent of this choice in the 
continuum limit. 

Case A . Here, the energy eigenvalue problem consists in: 

Hlfj^n) = -z(lW/27r)(>(0 n+1 ) - V#n)) = *Wn) , (54) 

with cf) n = 27m /TV, 1 < n < N, and the continuum limit will be considered momentarily. 
The complete orthonormal set of eigenfunctions and the eigenvalues are easily found: 

<M0n) = N- 1 ? 2 ex P [i(m + 5)<j) n ] , 1 < m < N , (55) 
E m = i(QN/2n){l-exp[2in{m + S)/N]) (56) 
Q{m + 5) , meN , (57) 

where 5 is an arbitrary real constant. Next, we turn to the symmetric discretization. 
Case B . Here we have instead: 

ftV#n) = -*(fiAY47r)(V(0 n+1 ) - V(0n-l)) = #V#n) , (58) 

with (f) n as before. This is solved by the same eigenfunctions as before, Eq. (|55jl . However, in 
this case the eigenvalues are real even before taking the continuum limit: 

E m = (£lN/2n)sm[2n(m + 5)/N] ^3° tt(m + 5) , m G N . (59) 

Thus, we learn that the spectrum in the continuum limit is indeed real and independent of the 
regularizations employed to suitably define the Hamiltonian, as it should be. 

The freedom to choose the constant S, which arises from the regularization of the Hamilton 
operator, is most wellcome. Choosing 5 = —1/2, we arrive at the quantum harmonic oscillator, 
having started from a corresponding classical system here. Thus, we recover in a straightforward 
way 't Hooft's result, derived from a cellular automaton See also Ref. [5] for the completion 
of a similar quantum model. In the following example we will encounter one more model of 
this kind and demonstrate its solution in detail. 

Also in the following example the regularized eigenvalues are complex and the real spectrum 
only emerges in the continuum limit. Again, this is due to the fact that we choose to discretize 
first-order derivatives asymmetrically, and could be avoided as shown. 

6.2 Quantum system with classical relativistic particle beneath 

Introducing proper time similarly as in Ref. j^j, the equations of motion and the Hamiltonian 
constraint of the reparametrization invariant kinematics of a classical relativistic particle of 
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mass m are given by: 

d T <f = m~V , d T p» = , (60) 
p ■ p — m 2 = , (61) 

respectively. Here we have ip a = (q°, . . . ,q 3 ;p°, . . . ,p 3 ), a = 1, . . . , 8; four-vector products 
involve the Minkowski metric, g^ v = diag(l, —1,-1,-1). 

Proceeding as before, we identify the effective Hamilton operator: 

U = -m~ l -n q ■ (p p = -m' l Tc q ■ (-id Wp ) , (62) 

corresponding to Eq. ()26|): the notation is as introduced after Eq. (|51|). however, involving four- 
vectors. Furthermore, the Hamiltonian constraint is represented by the operator: 

C = 5(0 p 2 - m 2 ) = 5{dl + m 2 ) , (63) 

according to Eqs. (J44j) and (J61|) . After a Fourier transformation, which replaces the variable ir q 
by a derivative (four-vector) +id x , and renaming ir p = x, we obtain: 

U = - m - 1 d x ■ d- x , C = 5(d? + m 2 ) , (64) 

i.e., a more transparent form of the Hamilton and constraint operators, respectively. 

Before embarking to further analyze this model, some general remarks seem in order here. 
It is well known from ordinary quantum (field) theory that the harmonic oscillator is peculiar in 
many respects. Therefore, the reader should not be misled by the results of Section 6.1, where we 
obtain a quantum harmonic oscillator spectrum from an underlying classical harmonic oscillator 
model. In particular, it may appear as if we have invented just one more quantization method, 
in line with quantization via canonical commutators, stochastic quantization, etc. However, we 
stress that this is not the case jHlE]- 

It seems an accident of the harmonic system that the usual quantized energy spectrum 
results here. This is revealed by the fact, already demonstrated in Refs. [31313, that localization 
with respect to the coordinate q of the underlying classical model has nothing to do with 
localization with respect to the operator q, which is introduced a posteriori when interpreting 
the emergent quantum Hamiltonian corresponding to said spectrum. Rather, such localized 
quantum (oscillator) states are widely spread over the q-space of the underlying classical model. 
We will encounter such operators in the following example again. Furthermore, the usual p, q- 
commutator algebra obtains corrections here, as long as the regularization is not removed. 

Therefore, generally, we do not expect to find the usual quantized counterpart of a classical 
reparametrization invariant model in the present approach based on discrete physical time. 
There will not be the usual one-to-one correspondence. To put it differently, the classical limit 
of emergent quantum theories cannot be expected to give back the underlying classical model. 
Further general remarks in this vein may be found in Ref. j^j. 

The following discussion of the free relativistic particle should be seen in this light. While 
any standard quantization method produces an unbound spectrum with notorious negative 
energy states, we show in detail here that careful application of the regularization indeed allows 
to produce an acceptable quantum model in the continuum limit, which is different from the 
underlying classical model. 
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The eigenvalue problem is solved after discretizing the system with a hypercubic (phase 
space) lattice of volume L 8 (lattice spacing I = L/N) and periodic boundary conditions, for 
example. Similarly as in the oscillator case, we here obtain the eigenfunctions: 

^k x ,kA x n, x n ) = N' 1 exp[i(k x + 8 X ) ■ x n + i(k x + 8 S ) ■ x n ] , (65) 

with coordinates x£ = ln^ and momenta k£ = 2irk' l /L, with 1 < n^, k^ < N, and where 5% are 
arbitrary real constants, for all \i = 0, . . . , 3 (analogously x%, fcf , 

These are the eigenfunctions, which are of plane-wave type, of the stationary Schrodinger 
equation, TiCip = ECip, with TC and C from Eqs. discretized analogously to Case A of 
Section 6.1. The four-vector momenta k£,k£ label different eigenfunctions, a familiar feature 
in quantum mechanics, and the constraint will be implemented shortly. Note that the phases 
8%, 8£ can still be chosen arbitrarily; however, since they must be fixed once for all, they cannot 
possibly absorb the variable four-momenta. Furthermore, we emphasize that the phase space 
lattice has been introduced for the only purpose of regularizing the Hamiltonian, to be followed 
by a suitable continuum limit. At present, however, it is unrelated to the underlying discreteness 
of physical time. 

The corresponding energy eigenvalues are: 

E hxM = -m-H- 2 ((exp[tl(k x + 5 x ) }-l)(exp[U(h + 5^°}-l) (66) 
- ]T(exp[zZ(^ + 8 x ) j ] - l)(exp[z/(^ + 8^] - 1)) 

= m -\k x + 5 x )-(k x + 5 s ) + 0(l) , (67) 

where is L is kept constant in the continuum limit, / — > 0; again, the four-momenta k£, k% 
label and determine the energies. Furthermore, in this limit, one finds that the Hamiltonian 
constraint requires timelike "on-shell" vectors k x , obeying (k x + 5 X ) 2 = m 2 , while leaving k x 
unconstrained. 

Continuing, we perform also the infinite volume limit, L — > oo, which results in a continuous 
energy spectrum in Eq. (|57jl. We observe that no matter how we choose the constants 8 X ,8 X , the 
spectrum will not be positive definite. Thus, the emergent model appears not to be acceptable, 
since it does not lead to a stable groundstate. 

However, let us proceed more carefully with the various limits involved and show that 
indeed a well-defined quantum model is obtained. For simplicity, considering (1 + 1) -dimensional 
Minkowski space and anticipating the massless limit, we rewrite Eq. (|6Tj) explicitly: 




where we suitably rescaled and renamed the constants and the momenta, which run in the range 
1 < k , k 0,1 < N = 2s + l. Furthermore, we incorporated the Hamiltonian (on-shell) constraint, 
such that only the positive root contributes: k° + 8° = {k 1 + 8 1 \ + 0(m 2 ) = — (k 1 + 8 1 ) +0(m 2 ). 
This can be achieved by suitably choosing 8 0,1 . 

In fact, just as in the previous harmonic oscillator case, the choice of the phase constants is 
crucial in determining the quantum model. Here we set: 
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We observe that this choice implies that in the continuum limit, with I — > 0, since iV — > oo, we 
have s — > oo and thus 6 1 — > — oo. Incorporating these phases, we obtain the manifestly positive 
definite spectrum: 



Efr, s /) = (^) 2 ((^ + s + \) + l) (>° + s + ~) + (si + s + \) + l) + 0( 



m) , (70) 



where also the (half)integer quantum numbers s z ^s z )l are introduced, all in the range — s < 
s z < s, which correctly replace fc 1 , A; ' 1 . 

Recalling the algebra of the ST/ (2) generators, with S z \s z ) = s z \s z ) in particular, we are led 
to consider the generic operator: 

h = S g + s + ~ , (71) 

i.e., diagonal with respect to |s 2 )-states of the (half)integer representations determined by s. 
In terms of such operators, we obtain the regularized Hamiltonian corresponding to Eq. (|7Uj) : 

H = (-^=-) 2 {l + h + h + ht + h(h + hi)) + 0(m) , (72) 



which will turn out to represent three coupled harmonic oscillators, including an additional 
contribution to the vacuum energy. 

A Hamiltonian of the type of h has been the starting point of 't Hooft's analysis j^j, which 
we adapt for our purposes in the following. 

Continuing with standard notation, we have S 2 = S 2 + S 2 + S 2 = s(s + 1), which suffices 
to obtain the following identity: 



'—i^ + Sy+l + h 2 ) " (73) 



h = 

2s + 

Furthermore, using S± = S x ± iS y , we define coordinate and conjugate momentum operators: 

q = l(aS- + a*S+) , p = ^(bS^ + b* S + ) , (74) 

where a and b are complex coefficients. Calculating the basic commutator with the help of 
[S + , S-] = 2S Z and using Eq. (fTTjl. we obtain: 

M = 1(1 ~ 2^Tl k) ' (75) 

provided we set Q(a*b) = — 2/ (2s + 1). Incorporating this, we calculate: 

S 2 + S 2 = (2S ± 1)2 (|a| 2 p 2 + \b\ 2 q 2 - (9fa • %b + Via ■ ^b){q,p}) . (76) 
In order to obtain a reasonable Hamiltonian in the continuum limit, we set: 

■ — , b=^L= , n^(^L-) 2 . (77) 
s + 1/2 Js + 1/2 V^L 



Then, the previous Eq. (|75j) becomes: 



nk = l? + l -n>e + w ± m (\v + m >) , (78) 
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reveiling a nonlinearly modified harmonic oscillator Hamiltonian, similarly as in Ref. j^j. 

Now it is safe to consider the continuum limit, 2s + 1 = N —>■ oo, keeping ^frnL and fl 
finite. This produces the usual g, ^commutator in Eq. (|75jl for states with limited energy and 
the standard harmonic oscillator Hamiltonian in Eq. ()78j) . 

Using these results in Eq. (172)1 . the Hamilton operator of the emergent quantum model is 
obtained: 

n = a + \ E (Pi + n% 2 ) + ^ + E (p/ + fi 2 g/) , (79) 

5=1,0,1 4ii J=0,1 

where the massless limit together with the infinite volume limit is carried out, m — > 0, L — > oo, 
in such a way that Q remains finite. 

The resulting Hamiltonian here is well defined in terms of continuous operators q and p, 
as usual, and has a positive definite spectrum. The coupling term might appear somewhat 
more familiar, if the oscillator algebra is realized in terms of bosonic creation and annihilation 
operators. 

In Refs. (HI El we calculated the matrix elements of operators q,p with respect to the SU(2) 
basis of primordial states in an analogous case, showing that localization of the quantum oscil- 
lator has little to do with localization in the classical model beneath, as we mentioned before. 

Finally, we remark that had we chosen 5 0,1 = 5 0,1 = 1/2, instead of Eqs. (|69|h then a relative 
sign between terms would remain, originating from the Minkowski metric, and this would yield 
the Hamiltonian 7i oc (1 + h)(h — hi), which is not positive definite. Similarly, any symmetric 
choice, <5 0,1 = 5 0,1 = 5 would suffer from this problem. 

This raises the important issue of the role of canonical transformations, and of symmetries 
in particular. It is conceivable that symmetries will play a role in restricting the apparent 
arbitrariness of the regularization. Independence from the choice of phase space coordinates 
employed in actual calculations might be a desirable feature. A preliminary study of an anhar- 
monic oscillator indicates that the simple discretization used here possibly needs to be improved, 
in order to fullfill this. 

6.3 Remarks 

Concluding this section, we may state that the features which we illustrated in the previous 
examples, especially the necessity of regularization (discretization), promise to make genuinely 
interacting models quite difficult to analyze. Interesting results may perhaps be found with the 
help of spectrum generating algebras and/or some to-be-developed perturbative methods. 

We find it interesting that our general Hamilton operator, Eq. ()26|) . does not allow for the 
direct addition of a constant energy term, while in its regularized form this is possible, due to 
the appearance of an arbitrary phase. Choosing the latter determines the groundstate energy, 
which survives the continuum limit. 

Leaving aside the Hamiltonian constraint momentarily, we observe that the Hamilton op- 
erator equation, also in cases with large numbers n of coupled degrees of freedom, amounts 
to systems of first order quasi-linear partial differential equations. They can be studied by 
the method of characteristics [21]. Thus, one finds one inhomogeneous equation, which can be 
trivially integrated. Furthermore, the remaining 2n equations for the characteristics present 
nothing but the classical Hamiltonian equations of motion. 

Therefore, integrable classical models can be decoupled at the level of the characteristic 
equations by canonical transformations, if they can be applied freely at the pre-quantum level. 
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Classical crystal-like models with harmonic forces, or free field theories, respectively, will thus 
give rise to corresponding free quantum mechanical systems. These are constructed in a different 
way in Refs. Presumably, the (fixing of a large class of) gauge transformations invoked there 
can be related to the existence of integrals of motion implied by integrability here. 

Finally, we point out that the Hamiltonian equations of motion preclude motion into classi- 
cally forbidden regions of the underlying system. Nevertheless, quantum mechanical tunneling 
is an intrinsic property of the quantum oscillator models that we obtained. 

7 Conclusions 

In the present paper, we pursue the view that quantum mechanics is an emergent description 
of nature, which possibly can be based on classical, pre-quantum concepts. 

Our approach is motivated by the ongoing construction of a reparametrization-invariant 
time. In turn, this is based on the observation that "time passes" when there is an observable 
change, which is localized with the observer. More precisely, necessary are incidents, i.e. ob- 
servable unit changes, which are recorded, and from which invariant quantities characterizing 
the change of the evolving system can be derived [HUH]- 

Presently, this has led us to assume the relation between the constructed physical time 
t and standard proper time r of the evolving system in the form of a statistical distribution, 
P(t; t) = P{r—t), cf. Eq. (0). We assume that the distribution is not explicitly time-dependent, 
which means, the physical clock is practically decoupled from the system under study. We 
explore the consequences of this situation for the description of the system. 

We have shown how to introduce "states", eventually building up a Hilbert space, in terms 
of certain functional integrals, Eqs. (jll|) - (fT2|) . which arise from the study of a classical gen- 
erating functional. The latter was introduced earlier in a different context, studying classical 
mechanics in functional form [221 I2H]- We employ this as a convenient tool, and modify it, in 
order to describe the observables of reparametrization-invariant systems with discrete time (Sec- 
tions). Studying the evolution of the states in general (Section 4), we are led to the Schrodinger 
equation, Eq. (J29|) . However, the Hamilton operator, Eq. (|26|) . has a non-standard form. 

We demonstrate that proper regularization of the continuum Hamilton operator is indis- 
pensable, in order to find a stable groundstate. Limitations imposed by symmetries or further 
constraints and consistency of the procedure need to be clarified. Other possible regularization 
schemes need to be explored. 

Coming back to the probabilistic relation between physical time and the evolution parameter 
figuring in the parameterized classical equations of motion: One would like to include the clock 
degrees of freedom consistently into the dynamics, in order to address the closed universe. This 
can be achieved by introducing suitable projectors into the generating functional [TH l2Uj. Their 
task is to represent quasi-local detectors which respond to a particle trajectory passing through 
in Yes/No fashion. In a more general setting, such a detector /projector has to be defined in 
terms of observables of the closed system. In this way, typical conditional probabilities can be 
handled, such as describing "What is the probability of observable X having a value in a range 
x to x + 5x, when observable Y has value y? v . Criteria for selecting the to-be-clock degrees 
of freedom are still unknown, other than simplicity. Most likely the resulting description of 
evolution and implicit notion of physical time will correspond to our distribution P(r; t) of 
Eq. (|7f). however, evolving with the system. Dissipative and memory effects, which will arise 
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after integration over clock degrees of freedom, might play a crucial role in deriving a unique 
large-scale quantum model. 

This stroboscopic quantization emerging from underlying classical dynamics certainly may 
be questioned in many respects. It might violate one or the other assumption of existing no-go 
theorems relating to hidden variables theories. However, we believe it is interesting to get 
closer to a working example, before discussing this. Unitary evolution and tunneling effects are 
recovered in this framework, under the proviso of regularization of the continuum formalism. 
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